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METHOD FOR ADAPTIVE IMAGE REGION PARTITION AND MORPHOLOGIC 
PROCESSING 



Technical Field 

This invention relates to a method for adaptive image region partition and morphologic 
5 processing. 

Background of the Invention 

Given an image that consists of a plurality of component regions, its image can be 
partitioned into zones of influence (ZOI). In the ZOI image, each pixel is associated with 
a component label. The component (foreground) pixels in the component regions are 
10 labeled by their component labels. The non-component (background) pixels are labeled 
by the labels of the component that is closest to the pixels. 

ZOI is a powerful tool to integrate results of multiple images acquired from different 
aspects of the same objects of interest. In one typical application scenario, the 

15 component regions are detected from one image or one detection mode. The component 
regions are used to form the ZOI. The ZOI is then used to condition the detection or 
perform measurements of other images or other detection modes. One application 
example is cell analysis in high content screening, cancer diagnosis, prenatal testing, 
genetic testing, tissue-based pathology assays, etc. In this application, cell nuclei are 

20 stained and imaged in one image for nuclear segmentation. Whole cells or specific cell 
characteristics are stained and acquired in another image for whole cell or specific cell 
characteristics segmentation. Since nuclei tend to be isolated from each other and whole 
cells tend to overlap, it is advantageous to divide the cell regions guided by the nuclear 
segmentation. ZOI could be used for cell region separation. 

25 

The ZOI operation is related to skeletal influence zones which are derived from Voronoi 
skeletons or medial axes that transform a 2D object into ID line representation (R. 
Ogniewicz and M. Hg, "Voronoi skeletons: Theory and applications," in IEEE Comp. 
Vision and Pattern Rec, June 1992, pp. 63-69; Ogniewicz, R. L. and Kubler. O: 
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30 Hierarchic Voronoi Skeletons, Pattern Recognition, nr. 28, pp. 343-359, 1995). 
Voronoi skeletons or medial axes are derived from the concept of distance 
transformation (A. Meijster, J. B. T. M. Roerdink, W. H. Hesselink, "A General 
Algorithm For Computing Distance Transforms In Linear Time", 2000; Mihail N. 
Kolountzakis, Kiriakos N. Kutulakos, "Fast Computation of the Euclidean Distance Map 

35 for Binary Images", Information Processing Letters, 1992). The ZOI operation is also 
related to recursive morphology (Haralick, R.M., Shapiro, L.G., "Computer and Robot 
Vision", Volume I, pp. 236-237, Addison-Wesley, 1992; Serra, J. "Image Analysis and 
Mathematical Morphology", PP. 385-387, Academic Press, 1982). The ZOI regions are 
constructed by recursive dilation and conditioning operations. 

40 

The distance transform can be computed sequentially by a two-pass procedure. The 
method is very fast. It requires 0(N 2 ) time for an N by N image. However, the 
skeletonization requires iterative thinning or similar techniques. The processing time of 
these methods is shape dependent and sometimes time consuming for complex shapes. 
45 The recursive morphology operation is not efficient for a general purpose computer and 
is time consuming for complex shapes. It is desirable to have a fast and time predictable 
method to partition regions into zones of influence. 

The prior art ZOI method uses the same distance metric for all components. This is 
50 useful only for simple applications where all The prior art ZOI method uses the same 
distance metric for all components. This is useful only for simple applications where all 
components are not differentiated. This is often not the case in practical applications. 
For example, in cell analysis applications, different cell types could be contained in the 
same image. In this case, it is disadvantageous to determine ZOI using the same distance 
55 metric across all cell types since large and small cells will be divided equally. It is 

desirable to have an adaptive ZOI algorithm that could adaptively apply a distance metric 
dependent on the component characteristics. For example, different distance scales could 
be used for different cell types or the distance could be scaled according to the size of the 
component. 



60 
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The same limitation exists in morphological dilation or erosion operations and their 
combinations, such as opening and closing. The prior art morphological operations use 
the same structuring element for all pixels in an image. The structuring element cannot 
be adaptively changed according to a pixel's component region property. Therefore, 
65 when trying to fill holes of a component region, other regions may be merged. Thus, the 
connectivity (homotopic relations) is destroyed. It is desirable to have an adaptive 
morphological processing method that adaptively adjusts the structuring element size and 
shape according to the characteristics of the components. 

70 Objects and Advantages 

This invention provides a method for adaptive image region partition and morphologic 
processing. The partition is computed sequentially by a two-pass procedure similar to 
distance transformation. In addition to the shortest distance, the method tracks, 
propagates, and updates the component label where a pixel's shortest distance is derived. 

75 The method is very fast and the processing time is predictable and is component shape 
independent. It requires 0(N 2 ) time for an N by N image and is ideal for a use in a 
general purpose CPU. It could adaptively apply a distance metric depending on 
component characteristics. It enables adaptive morphological processing that allows 
adaptive adjustment of the structuring element size and shape according to characteristics 

80 of the components. Furthermore, the method enables adaptive integration of multiple 
channel information in cell analysis. 

The primary objective of the invention is to perform fast image region partitioning. A 
secondary objective is to achieve processing speed mostly independent of component 

85 shape. Another objective of the invention is to adaptively apply a distance metric 
dependent on component characteristics. The fourth objective is to provide adaptive 
morphologic processing that allows adaptive adjustment of the structuring element size 
and shape according to the characteristics of the components. The fifth objective of the 
invention is to provide adaptive integration of multiple channel information in cell 

90 analysis for improved image measurement accuracy. The sixth objective of the invention 
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is to reduce cost and increase availability by facilitating real time applications using 
general purpose CPUs. 



Summary of the Invention 

A fast image region partition method receives a component labeled image and performs a 
95 two pass Zone Of Influence (ZOI) creation method to create a Zone Of Influence (ZOI) 
image. The two pass ZOI creation method performs a first pass scan to create a first pass 
intermediate distance image and a shortest distance component label image. It then 
performs a second pass scan using the first pass intermediate distance image and the 
shortest distance component label image to create a background distance transform image 
100 and a updated shortest distance component label image. An adaptive image region 

partition method receives a component labeled image and performs an adaptive two pass 
ZOI creation method to create an adaptive ZOI image. The distance lengths of the two 
pass adaptive ZOI creation method depend on their associated component labels. 

105 An adaptive cell segmentation method receives a nuclei mask image and a cell mask 

image. It performs adaptive nuclei region partition using the nuclei mask image to create 
adaptive nuclei mask ZOI. An adaptive cell region separation method uses the cell masks 
and the adaptive nuclei mask ZOI to generate adaptive cell separated regions. 

1 10 An adaptive dilation method receives an image and performs an adaptive background 
distance transform to create an adaptive background distance transform image. A 
threshold is applied to the adaptive background distance transform image to generate 
adaptive dilation image output. An adaptive erosion method receives an image and 
performs an adaptive foreground distance transform to create an adaptive foreground 

115 distance transform image. A threshold is applied to the adaptive foreground distance 
transform image to generate adaptive erosion image output. 
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Brief Description of the Drawings 

120 The preferred embodiment and other aspects of the invention will become apparent from 

the following detailed description of the invention when read in conjunction with the 

accompanying drawings, which are provided for the purpose of describing embodiments 

of the invention and not for limiting same, in which: 

Figure 1 A shows a Euclidean disc; 
125 Figure IB shows the Euclidean distance transformed image of 1A; 

Figure 1C shows the city block distance transformed image of 1A; 

Figure ID shows the chessboard distance transformed image of 1 A; 

Figure 2 shows the medial axis of a rectangle defined in terms of bi-tangent circles; 

Figure 3 shows the processing flow for the fast zones of interest creation method; 
130 Figure 4 shows the processing flow for a cell segmentation method; 

Figure 5A shows an illustrative image containing 5 cells of different types; 

Figure 5B shows the nuclei masks of 5A; 

Figure 5C shows the cell masks of 5A; 

Figure 5D shows the overlay of nuclei mask ZOI boundaries on the cell image of 5 A; 
135 Figure 5E shows the distance transform of the cell masks of 5C; 

Figure 5F shows the overlay of adaptive nuclei mask ZOI boundaries on the cell image of 
5A; 

Figure 5G shows the cell region segmentation of 5A after removing the pixels along the 
adaptive nuclei mask ZOI boundaries. 

140 Detailed Description of the Invention 

I. Distance Transform 

Given an image I, the foreground distance transform assigns each foreground pixel the 
shortest distance between the pixel and a background pixel. Conversely, the background 
145 distance transform assigns each background pixel the shortest distance between the pixel 
and a foreground pixel. 
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The distance transform is computed sequentially by a two-pass method. The first 
(forward) pass scans in a left-right top-bottom raster scan order. The second (backward) 
150 pass scans in a reverse right-left bottom-top order. In the first pass for a foreground 

distance transform, the output U(x,y) at pixel position (x,y) is determined by the value of 
input image I and the values at previously computed positions of U by 



155 



U(x,y) = { 



min {U (i f j) + /(/, ;)}; V/(x, y) e foreground 



0\ V/(x,/y) e background 



In one embodiment, the boundary condition is handled by setting U(x,-1) = 0 and U(-l,y) 
= 0. N B (x,y) is the backward neighbors of the pixel (x,y). The backward neighbor 
contains a selected set of adjacent neighbors that are already scanned in the first pass. 
l(i,j) is the neighbor position dependent distance length. When l(i,j) = 1 for all (i, j) and 
160 N B (x,y) is the upper half of the 4-connected neighbor, thus, N B (x,y) = {(x-l,y), (x,y-l)}, 
then the resulting distance is the well-known 4-connected city block distance (Li). When 
l(i,j) = 1 for all (i, j) and N B (x,y) is the upper half of the 8-connected neighbor, thus, 
N B (x,y) = {(x-l,y), (x-l,y-l), (x,y-l), (x+l,y-l)},then the resulting distance is the well- 
known 8-connected chessboard distance (Loo). For the 8-connected neighbor, when l(x- 

165 l,y)=l, l(x-l,y-l)= V2 , l(x,y-l)=l, and l(x+l,y-l)= V2 , the resulting distance is the 
well-known Euclidean distance (L 2 ). 

The second (backward) pass creates the distance transform image D by 
170 D(x, y) = min{ min (D(i, l(U j)),U (*, y) } 

(i,j)eN F (x,y) 

N F (x,y) is the forward neighbors of the pixel (x,y). The forward neighbor contains a 
selected set of adjacent neighbors that are already scanned in the second pass. For 4- 
connected city block distance, N F (x,y) is the lower half of the 4-connected neighbor. 
175 That is, N F (x,y) = {(x+l,y), (x,y+l)}. For 8-connected chessboard distance and 
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Euclidean distance, N F (x,y) is the lower half of the 8-connected neighbor. That is, 
N F (x,y) = {(x+l,y). (x+l,y+l), (x,y+l), (x-l,y+l)}. 



The images in Figure 1 A to Figure ID illustrate the three distance transforms applied to a 
Euclidean disc. Figure 1 A shows the Euclidean disc; Figure IB shows the Euclidean 
180 distance transformed image of Figure 1 A; Figure 1C shows the city block distance 

transformed image of Figure 1A; Figure ID shows the chessboard distance transformed 
image of Figure 1A. 

II. Medial Axis 

185 Medial axis or skeleton is initially introduced for shape description. The following 

hypothetical scenario can help visualize a medial axis. Let the shape be a region of grass 
surrounded by bare soil. Simultaneously at all boundary points of the grass region light a 
fire that burns uniformly from the border of the grass region to its interior. Eventually 
the firs burning from one part of the border will meet the fire burning from another part 

190 of the border. Since all the grass will have been burned at the meeting place, the fire will 
quench and be extinguished. The locus of points on the arc of fire extinction is the 
medial axis of the region. Another way to think about the medial axis is as the loci of 
centers of bi-tangent circles that fit entirely within the foreground region being 
considered. Figure 2 illustrates this concept for a rectangular shape. The medial axis of a 

195 rectangle 200 is defined in terms of bi-tangent circles in Figure 2. 

Each point on the medial axis has an associated value that is proportional to the time it 
took for the fire to reach the given point from the time the grass fire was set. The medial 
axis with its medial axis distance function is called the medial axis transform. It is an 
200 information preserving representation of shape. To see this, just consider running the 
grass fire backward. With time running in reverse, set a grass fire on each point of the 
medial axis exactly at the time the original grass fire is extinguished at that point. The 
boundary of the fire at time t=0 would be the boundary of the original given shape. 



- 8/21 - 

205 The medial axis can be produced in two main ways. The first is to use morphological 
thinning that successively erodes away pixels from the boundary (while preserving the 
end points of line segments) until no more thinning is possible, at which point what is left 
approximates the medial axis (skeleton). An alternative method is to first calculate the 
distance transform of the image. The skeleton then lies along the singularities (i.e. creases 

210 or curvature discontinuities) in the distance transform image. The medial axis is often 
described as being the 'locus of local maxima' on the distance transform image. If the 
distance transform image is displayed as a 3-D surface plot with the third dimension 
representing the grayscale value, the medial axis can be imagined as the ridges on the 3-D 
surface. 

215 

m. ZONES OF INFLUENCE 

The distance transform can be applied to the background of the image and the skeleton or 
medial axis of the background could be defined. One type of background skeleton is the 
220 zone of influence. To define the zones of influence, begin with a representation of a set S 
as the union of its components Cj, . . ., Cn: 

N 

s = Uc, 

225 The zone of influence 7^ associated with component C n is the set of all points that are 
closer to C n than to any other components. That is, 

Z„ = {P\ d(p,C n )< d(p,CJ,m * n) 

230 where d(p, C) denotes the distance from the point p to the set C. 

Given an image consisting of a plurality of component regions, its background (non- 
component regions) can be partitioned into zones of influence (ZOI). In the ZOI image, 
each pixel is associated with a component label. The foreground pixels in the component 
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235 regions are labeled by their component labels. The background pixels are labeled by the 
labels of the component that is closest to the pixel. Note that the zones of influence could 
be defined even though the pixels in a component are not connected. 

We now describe the fast image region partition (zones of influence creation) method of 
240 this invention. The processing flow of the fast zones of interest creation method is shown 
in Figure 3. As shown in Figure 3, the input image 300 is subject to component labeling 
306 that assigns a unique label to each component; this results in the component labeled 
image 302. The component labeled image 302 is used by the two pass ZOI creation 308 
method that uses a distance transformation-like two pass method to create the ZOI image 
245 304. The method is very fast. It requires 0(N 2 ) time for an N by N image. 

EI. 1 Component labeling 

In the case that the connectivity is required for each component, the standard connected 
250 component labeling method could be used to generate the component labeled image. As 
mentioned previously, pixels in a component do not have to be connected. Those skilled 
in the art should recognize that other labeling methods such as manual labeling, semi- 
automatic labeling, image segmentation and clustering could be used. 

255 m.2 Two Pass ZOI Creation 

Given a component labeled image L, the two pass ZOI creation method sequentially 
computes the ZOI image by a two-pass method. Similar to the distance transformation 
method, the first (forward) pass scans in a left-right top-bottom raster scan order. The 
260 second (backward) pass scans in a reverse right-left bottom-top order. In the first pass, 
an output first pass intermediate distance image value F(x,y) and a shortest distance 
component label image Z(x,y) at pixel position (x,y) are determined by the value of 
component labeled image L and the previously computed F and Z values by the following 
algorithm: 



265 



F(x,y) = < 



Z(x,y) = 
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min F(i, j) + j); \/L(x, y) e background 

(i,j)sN B (x,y) 



0;VL(;t, y) e foreground 



Z(arg min {F(i, j) + j)}); VL(jc, y) e background 

(iJ)eN B (x,y) 



L(x, y); \/L(x, y) e foreground 



270 The first pass intermediate distance image F stores the first pass intermediate distance 
values and the shortest distance component label image Z(x,y) stores the label of the 
component that has the shortest intermediate distance to pixel (x,y). 

In one embodiment, we handle the boundary condition by setting F(x,-1) = 0 and F(-l,y) 
275 = 0. Those skilled in the art should recognize that other methods of boundary condition 
handling could be used and they are all within the scope of this invention. 

The second pass produces the background distance transform image D and updates the 
shortest distance component label image, Z(x,y), by 

280 

D(x, y) = min{ min (D(i, j) + j)), F(x, y)} 

OJ)eN F (x,y) 



Z(x,y) = < 



Z(arg min {D(i,j) + l(i,j)));If min {D(i,j) + l(i,j)}<F(x,y) 

VJ)eN h (x,y) (i,j)eN F (x,y) 

Z(x, y);Otherwise 



285 



The updated Z(x,y) of the second pass stores the label of the component that has the 
shortest distance to pixel (x,y). Therefore, the final Z is the resulting ZOI image. The 
length function l(i,j) and the forward and backward neighbors, N F (x,y) and N B (x,y), of a 
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pixel (x,y) can be selected so that the distance metric used for the ZOI could be city block 
290 distance (Li), chessboard distance (Loo), Euclidean distance (L 2 ), or other distance 
metrics. 

Those skilled in the art will recognize that the scanning directions of the first pass and the 
second pass could be reversed. Furthermore, other complementary scanning directions 
295 could be used and they are all within the scope of this invention. For example, the first 
pass could be scanned in a right-left top-bottom scan and the second pass could be 
scanned in a left-right bottom-top scan order. Other configurations of the forward and 
backward neighbors can also be used. 

300 IV. Adaptive Image Region Partition 

The prior art zone of influence method assumes that a fixed distance metric is used 
throughout the ZOI determination. This treats each component equally and represents a 
major limitation in practice. For example, in cell analysis applications, different cell 

305 types could be contained in the same image. In this case, it is disadvantageous to 

determine ZOI using the same distance metric across all cell types since large and small 
cells will be divided equally. It is desirable to have an adaptive ZOI algorithm that could 
adaptively apply a distance metric depending on component characteristics. For 
example, different distance scales could be used for different cell types or the distance 

310 could be scaled according to the size of the component. 

The adaptive image region partition method of this invention is accomplished by making 
the length function l(i,j) dependent on component labels. The first pass of the adaptive 
two pass ZOI creation method outputs a first pass intermediate adaptive distance image 
315 F A and an adaptive shortest distance component label image Z A by the following 
algorithm 
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min { F A (i, j) + , j), Z A (i, j)] } ; \/L(x, y) e background 

(iJ)eN B (x, y ) 



0;\/L(x,y) e foreground 



320 



Z A (x,y) = < 



Z A (arg min {F A (i, j) + j),Z A (i, j)]}); VL(*, y) e background 



(i,j)eN B (x,y) 



L(x, y);VL(x, y) s foreground 



325 



The second pass produces the adaptive distance transform image D A and updates the 
adaptive shortest distance component label image, Z A (x,y), by 

D A (x, y) = min { min (D A (i, j) + /[(/, j), Z A (i, j)]), F A (x,y)} 

(iJ)eN (x,y) 



Z A (x,y) = ' 



Z^arg min {D A (i,j) + l[(i,j),Z A (i,j)]}),If min {D A (i, j) + j),Z A (i, j)]} < F A (x,y) 
Z A (x, y) y Otherwise 



330 



335 



Note that the length l[(i,j),Z A (ij)] is a function of the neighboring position (i j) as well as 
the neighbor's associated component labels, ZaOj). When the forward and backward 
neighbors, N F (x,y) and N B (x,y), are defined as a super set of all desirable neighbors, a 
desirable distance metric for a component label 1 could be constructed by providing 
appropriate length values among the neighbor points. If the super set neighbor includes 
unnecessary positions, they can be effective disabled by the assignment of large length 
values for the positions. 



340 



Those skilled in the art will recognize that the scanning directions of the first pass and the 
second pass could be reversed. Furthermore, other complementary scanning directions 
could be used and they are all within the scope of this invention. For example, the first 
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pass could be scanned in a right-left top-bottom scan and the second pass could be 
scanned in a left-right bottom-top scan order. Other configurations of the forward and 
backward neighbors can also be used. 

345 

In one embodiment of the invention, a component specific weighting factor, W(c), is 
applied to the original length function for adaptive image region partition. That is, 

l[(i,j),Z(i,j)] = l(i,j)*w(Z(i,j)) 

350 

Where w(c) could be related to the size of the component or other factors such as the 
confidence or the types of the components. 

Those skilled in the art should recognize that the length function could be any linear or 
355 nonlinear functions of the component labels. They are all within the scope of this 
invention. 

The image region partition method can be used for cell analysis. Figure 4 shows the 
processing flow of a cell segmentation method. As shown in Figure 4, the input nuclei 
360 mask image 400 is processed by the nuclei region partition 408 method to create nuclei 
mask ZOI 404. The nuclei region partition 408 can be performed using the fast zones of 
interest creation method. The cell mask image 402 and the nuclei mask ZOI 404 are used 
by a cell region separation 410 step to create cell separated regions 406 output. 

365 The cell region separation 410 step considers the cell mask having the same component 
label in the nuclei mask ZOI 404 as one cell region. This results in the cell separated 
regions 406 output. 

Figure 5A to Figure 5D illustrates an example of cell segmentation method. Figure 5 A 
370 shows an illustration of an image with 5 cells of different types. The component image 
contains nuclei masks and the whole cell image contains cell masks as illustrated in 
Figure 5B and Figure 5C. 
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Figure 5D illustrates the overlay of the nuclei mask ZOI using static Euclidean distance. 
375 The cell mask can be separated by the nuclei mask ZOI as shown in Figure 5D. Note that 
the static distance cuts the cell in the middle of two nuclei, which incorrectly partitions 
part of cytoplasm of the large cells into smaller cells. 

The adaptive image region partition method could be applied to improve performance. 
380 An adaptive cell segmentation method inputs a nuclei mask image and a cell mask image. 
It performs adaptive nuclei region partition using the nuclei masks to create adaptive 
nuclei mask ZOI. It then performs adaptive cell region separation using the cell masks 
and the adaptive nuclei mask ZOI to generate adaptive cell separated regions output. 

385 The adaptive nuclei region partition method uses cell size estimate for the weighting 

factor of the length function. Figure 5E shows the distance transform of the cell masks. 
In one embodiment of the invention, the cell size estimate is performed by averaging the 
cell distance values within a nucleus. Let the cell size estimate for nucleus i be Si. In 
one embodiment of the invention, the median value of Si among all nuclei in the image is 

390 determined as S'. A weighting factor could be determined for each component as w(i) = 
Min (a, Max( b, S VSi)) Where a > b. a and b are determined by the expected size 
difference between cells in the same image. In one embodiment of the invention, the 
values are set to a = 2.0 and b = 0.5. Figure 5F shows the overlay of the adaptive nuclei 
mask ZOI boundaries on the cell image. An enhancement could be performed by 

395 excluding the pixels close to the adaptive nuclei mask ZOI boundaries. This further 

removes incorrect cytoplasm pixels to allow accurate cytoplasm measurements. Figure 
5G shows the results of this exclusion. 

V. Adaptive Morphological Processing 

400 

The morphological dilation and erosion operations can be accomplished by distance 
transform followed by thresholding. Dilation can be accomplished by applying a 
distance transform on the image background and then applying the thresholding. The 
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thresholding operation selected the pixels whose values are less than or equal to the 
405 threshold value in the resulting mask. Erosion can be accomplished by applying a 

distance transform on the image foreground and then applying the thresholding. In this 
method, the distance metric determines the structuring element shape and the threshold 
value determines the size of the structuring element. 

410 The prior art morphological dilation, erosion operations and their combinations such as 
opening and closing use a common structuring element for all pixels in an image. The 
structuring element cannot be adaptively changed according to a pixel's component 
region property. Therefore, when trying to fill the holes of a component region, other 
regions may be merged and therefore, the connectivity (homotopic relations) destroyed. 

415 

As described in section IV, the adaptive image region partition of this invention generates 
not only the Zone of influence image Z but also the adaptive background distance 
transform image D A . Applying threshold to the adaptive background distance transform 
image results in the adaptive dilation image output. 

420 

Adaptive erosion can be accomplished by performing adaptive foreground distance 
transform to create an adaptive foreground distance transform image. This is followed by 
applying a threshold on the adaptive foreground distance transform image to create the 
adaptive erosion image output. 

425 

In one embodiment of the invention, the adaptive foreground distance transform is 
accomplished using a two-pass procedure. Similar to the standard distance 
transformation method, the first (forward) pass scans in a left-right top-bottom raster scan 
order. The second (backward) pass scans in a reverse right-left bottom-top order. In the 
430 first pass, the output image values F(x,y) at pixel position (x,y) is determined by the value 
of component labeled image L and the values at previously computed positions of F by 
the following algorithm: 



435 



F(x,y) = < 
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min {F(iJ) + l[(iJ) 9 L(iJ)]} m yL{x 9 y) e foreground 
0; VL(x, g background 
The second pass produces the adaptive foreground distance transform image D F (x,y), by 



D F O, y) = min{ min (D F (i, j) + /[(/, L(i, ;)]), F(x, y) } 

(i,j)eN F (x,y) 



440 Note that the length l[(i,j),Z(i,j)] is a function of the neighboring position (i,j) as well as 
the component labels, L(ij). When the forward and backward neighbors, N F (x,y) and 
N B (x,y), are defined as a super set of all desirable neighbors, a desirable distance metric 
for a component label 1 could be constructed by providing appropriate length values 
among the neighbor points. If the super set neighbor includes unnecessary positions, they 

445 can be effective disabled by the assignment of large length values on these positions. 

Those skilled in the art will recognize that the scanning directions of the first pass and the 
second pass could be reversed. Furthermore, other complementary scanning directions 
could be used and they are all within the scope of this invention. For example, the first 
450 pass could be scanned in a right-left top-bottom scan and the second pass could be 

scanned in a left-right bottom-top scan order. Other configurations of the forward and 
backward neighbors can also be used. 

The invention has been described herein in considerable detail in order to comply with 
455 the Patent Statutes and to provide those skilled in the art with the information needed to 
apply the novel principles and to construct and use such specialized components as are 
required. However, it is to be understood that the inventions can be carried out by 
specifically different equipment and devices, and that various modifications, both as to 
the equipment details and operating procedures, can be accomplished without departing 
460 from the scope of the invention itself. 



